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ABSTRACT 

Performing ground-based submillimetre observations is a difficult task as the mea¬ 
surements are subject to absorption and emission from water vapour in the Earth’s 
atmosphere and time variation in weather and instrument stability. Removing these 
features and other artifacts from the data is a vital process which affects the charac¬ 
teristics of the recovered astronomical structure we seek to study. In this paper, we 
explore two data reduction methods for data taken with the Submillimetre Common- 
User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). 
The JCMT Legacy Reduction 1 (JCMT LRl) and The Gould Belt Legacy Survey 
Legacy Release 1 (CBS LRl) reduction both use the same software (Starlink) but dif¬ 
fer in their choice of data reduction parameters. We find that the JCMT LRl reduction 
is suitable for determining whether or not compact emission is present in a given re¬ 
gion and the CBS LRl reduction is tuned in a robust way to uncover more extended 
emission, which better serves more in-depth physical analyses of star-forming regions. 
Using the CBS LRl method, we find that compact sources are recovered well, even at 
a peak brightness of only 3 times the noise, whereas the reconstruction of larger ob¬ 
jects requires much care when drawing boundaries around the expected astronomical 
signal in the data reduction process. Incorrect boundaries can lead to false structure 
identification or it can cause structure to be missed. In the JCMT LRl reduction, 
the extent of the true structure of objects larger than a point source is never fully 
recovered. 

Key words: ISM: structure - techniques: image processing - stars: formation - 
submillimetre: ISM 


1 INTRODUCTION 


In an effort to probe the cold dust in several well-known 
nearby star-forming regions with the James Clerk Maxwell 
Telescope (JCMT), the Gould Belt Survey (CBS, |Ward-| 
Thompson et al.|2007bl has been performing submillimetre 


continuum observations using the Submillimetre Common- 
User Bolometer Array 2 (SCUBA-2). The SCUBA-2 instru¬ 


ment (Holland et al. 20131 is a wide field submillimetre 


bolometer camera with an unprecedented 10,000 pixels. The 
device maps regions at 450 pm and 850 pm simultaneously 
with effective resolutions of 9.6” and 14.1" (Dempsey et al. 


20131, respectively. Mapping the sky over 100 times faster 
than its predecessor, SCUBA ( [Holland et al.||1999p , this next 
generation detector has offered star formation researchers 
the chance to analyse nascent stellar systems at submil¬ 
limetre wavelengths with a higher efficiency and a broader 
context than ever before with a single dish, ground-based 
telescope (see, for example, Sadavoy et al.|[2013 Salji et al. 


|2015a| and jPattle et al.|2015p . 

Reducing SCUBA-2 data is a complex process with a 
variety of solutions, each designed to best uncover particu¬ 
lar features (e.g. bright, compact emission versus faint, dif¬ 
fuse emission). This variety is, in part, due to the nature of 
the observations. While performing observations using the 
PONG scanning mode (the method employed in this study; 
see jKackley et al.|2010| and |Holland et al.|2013| ), the JCMT 
continuously scans back and forth across the sky at differing 
angles to HU in a circular pattern. By visiting the same lo¬ 
cations from different angles and at different times, sources 
of low frequency drift are manifested as a length-scale fea¬ 
ture which can be separated from the sources of interest so 
that spatially invariant structures can be identified. These 


et al.|20T3 1 but are dominated by atmospheric noise, which 
varies temporally. In this way, many individual bolometers 
will observe each sky location. Constructing a final image re¬ 
quires identifying significant structure and removing those 
large-scale features created by the noise. Quantifying how 
the observing process affects the observed signal is therefore 
of utmost importance as it characterises how confident we 
can be that we are detecting true astronomical emission and 
not artificial constructs at each scale of interest. 

Determining the optimal image reconstruction of sub¬ 
millimetre bolometer data has been an area of acute inter¬ 
est for researchers using JCMT data since before the op¬ 
eration of the original SCUBA instrument. For example. 


Richer (19921 presented the “maximum entropy” reconstruc¬ 


tion method, which is based on the assumption that the 
most likely reconstruction is the one which has maximum 
entropy (see Narayan fc Nityananda|1986 for a review). Wil 


[son et ah] ( |1999| ) an d]P ierce-Pricej ( |200^ went on to develop 
maximum entropy reconstruction methods for data taken 
by SCUBA. Several other methods were also employed. For 
instance, |Borys| (|2002|) wrote a SCUBA map-making algo¬ 
rithm based on a least squares approach while [Johnston^ 


et al. (2000a I developed a “matrix inversion” reconstruction 


which accounted for the varying uncertainty in a given map’s 
background noise. The latter was based on a technique for 
the Wilkinson Microwave Anisotropy Probe (WMAP) de- 
scrib ed by [Wright et al.|(|1996p . At approximately the same 
time, Jenness et al. (20001 implemented a Fourier deconvo¬ 


lution method developed by Emerson (19951 (Emerson 2) as 


sources of drift include instrument-based noise (see Holland 


the reduction method for the standard SCUBA data analysis 
software at the JCMT. 

The issue of effective and accurate submillimetre data 
reduction is not unique to the JCMT. For example, the Bolo- 
cam Galactic Plane Survey (BGPS; see Aguirre et al.|[2MT 
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for the 1.1mm survey and [Merello et aL]|2015| for the 350 
/rm images) uses an algorithm called CRUSH (The Com¬ 
prehensive Reduction Utility for SHARC-II; Kovacs|2006| ) to 
reduce the 350 /rm continuum images taken using the Sub¬ 
millimetre High Angular Resolution Camera II (SHARC-II; 
Dowell et al. 20031 at the Caltech Submillimetre Observa¬ 


tory (CSO). Like SCUBA-2, SHARC-II is a ground-based 
bolometer array which must contend with instrumental im¬ 
perfections, atmospheric interference, and electronic noise. 
The CRUSH algorithm is an iterative process which deter¬ 
mines a series of gain and weight estimations for each array 
component one-by-one while filtering noise spikes and other 
bad data (such as cosmic rays and electronic discharges; 
Kovacs|[2006 for more information). Other examples in¬ 


clude SANEPIC, the maximum-likelihood map-making al¬ 
gorithm used for the Balloon-borne Large Aperture Submil¬ 


limeter Telescope (BLAST; Patanchon et al. 20081 experi¬ 
ment which uses a series of approximations to reduce the re¬ 
quired computational power necessary for image production, 
and the AzTEC (Aztronomical Thermal Emission Camera, 
located at the Large Millimetre Telescope) Data Reduction 
Pipeline, which uses Principal Component Analysis (PCA) 
to model and remove correlated components of the bolome¬ 
ter signals. The latter method, however, is only ideal for the 
recovery of compact structures as extended emission will 
have correlated components that will be flagged as noise 
I Scott et al.||2008 K 

In this paper, we explore the data reduction method 
used by the JCMT CBS team for their first Legacy Release 
(CBS LRl; Kirk et ah, in prep.) of SCUBA-2 data by first 
comparing and contrasting with another method employed 
by the JCMT science and computing group (JCMT LRl, 
Graves et al., in prep.) then analysing how well each method 
preserves structures with known properties. We focus here 
only on the 850 fim maps. The same arguments can be ex¬ 
tended to 450 /rm maps but there is less high signal-to-noise 
data available for comparison and the absolute calibration 
uncertainties are higher. 

To construct an image from the raw SCUBA-2 data, 
the CBS LRl and the JCMT LRl reduction teams both use 
the makemap algorithm found in Starlink’i0 smurf pack¬ 
age (Chapin et al. 2013 also see the SCUBA-2 data re- 
duction cookbook [Thomas fc Currie|[20T^ . The makemap 
algorithm employs an Iterative technique which applies flat 
field corrections, performs an extinction correction, models 
and removes noise correlated between detectors, estimates 
and masks the astronomical signal, and measures the noise 
of each bolometer contributing to each pixel in the reduced 
image. In total, for all the different observational strategies 
available, makemap has over one hundred configurable pa¬ 
rameters which are supplied in a text file called the “Dy¬ 
namic Iterative Map Maker configuration file”, or, “dimm- 
config file”. With these parameters, the user has significant 
control over the entire reduction algorithm from the pre¬ 
processing stages to what is included in each model (as¬ 
tronomical, noise, etc.). Varying the parameters supplied in 
makemap's dimmconfig file, therefore, will give rise to dif¬ 
ferent final images. Many of these differences will be subtle 


http://www.starlink.ac.uk 


but any differences between reductions should be well un¬ 
derstood. 

The dimmconfig file employed Is dependent on the sci¬ 
entific goals desired. For example, since the beginning of the 
SCUBA-2 mapping initiative, the JCMT Science Archive’s 
( [Economou et al.|2015[ ) Legacy Release project has been fo¬ 
cused on producing public images and subsequent catalogues 
which could answer the simple question: “Is there compact 
emission in an observed region?”. In this way, the archive can 
provide a useful set of information to a wide variety of user 
projects without doing more complex analyses that require 
careful human oversight. 

To this end, the JCMT LRl data reduction was tuned 
to downgrade extended structure while effectively identify¬ 
ing compact emission. The CBS SCUBA-2 team, however, 
has a wide range of specific scientific interests including but 
not limited to an investigation of the emissivity spectral in¬ 
dex of the thermal dust emission (Hatchell et al. 2013| Coude 
et al., in prep.) and its relationship to temperature and col¬ 
umn density (in conjunction with Herschel Space Telescope 
data; Sadavoy et al. 2013 Chen et ah, in prep.), structure 
mapping (Salji et al. 2015a Salji et al. 2015b[), in vestiga- 
tions into protoplanetary disks ( [Dodds et al.| 2015[ |Buckle| 
et al. 2015| Broekhoven-Fiene et ah, in prep.), fragmenta¬ 


tion analyses (Mairs et al., in prep.), virial studies (Rattle 
|et al. [^15|), clustering (Lane et al., in prep., Kirk et al., in 


prep.), and radiative feedback processes around young stars 
(Rumble et al. 20151. The CBS LRl is the most current 


data reduction product available from the CBS team to suit 
these individual goals while remaining consistent across all 
the star-forming regions observed by the CBS. This broad 
mandate means that the data reduction strategy requires 
both compact and extended emission recovery. The latter 
Is a more difficult process as the separation of large, faint 
structures from time-varying noise is non-trivial. To this 
end, the CBS LRl reduction employs a two-step process 
in which significant recovered emission is first identified au¬ 
tomatically based on the signal-to-noise ratio (SNR) then, 
based on those resulting maps, a rigorous analysis is per¬ 
formed wherein the observer uses the structure detected and 
in some cases knowledge from other data sets (for example, 
from the Herschel Space Observatory) to define a boundary 
around any emission that is likely real. After the significant 
structure is more rigorously identified by the observer, the 
map undergoes a second round of data reduction. 

The format of this paper is as follows: In Sectionwe 
outline in more detail the specific differences between the 
JCMT LRl reduction and CBS LRl reduction methods. In 
Section we compare the final maps produced using each 
method in three regions of interest within the Orion A South 
star-forming complex. In Section we explore the scales 
and flux levels at which the reductions preserve structure 
by constructing artificial Gaussian sources and recovering 
them after they have been processed by makemap in a pure 
noise field. In Sectionwe investigate the effect of changing 
the flux density threshold at which SCUBA-2 data are con¬ 
sidered to be significant astronomical signal. We then carry 
on this investigation by determining the effects of changing 
the size of a user-defined boundary around emission that 
has been labelled as significant. In Sectionwe give a brief 
overview of how data reduction can affect the results of com¬ 
mon observational metrics used when studying star forming 
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regions. Finally, in Section]^ we summarise the main results 
and present our conclusions. 


2 DATA REDUCTION PARAMETERS 


The SCUBA-2 data reduction process, makemap, is ex¬ 
plained in detail by Chapin et al. (20131 (see their Fig¬ 
ure 6). To summarise here, the raw data from each scan 
is hrst assembled into a continuous time series and a flat 
held correction is applied. Then, the data is downsampled 
and discontinuities such as steps, spikes, and gaps are re¬ 
paired. Following that, the mean of each bolometer time 
series is removed from the map and the iterative portion of 
the mapmaking procedure begins. 

There are six models iteratively constructed and Im¬ 
proved upon in hve steps before the hnal map is created: 


1. The COM and GAI models hrst estimate and remove 
the common mode signal across the bolometers at each 
time step, respectively. The common mode removal can be 
performed over the full SCUBA-2 focal plane array (~400” 
scale) or over each sub-array individually (~200” scale; 
see the com.perarray parameter in Sections |2.1| and |2.2[ 
below). 


2. The EXT model then corrects for extinction. 


3. Next, the FLT model (based on a high-pass hltering al¬ 
gorithm) removes independent low-frequency noise associ¬ 
ated with each individual bolometer directly from the time 
series. The physical scale to which we hlter is initially input 
as a length (see the f It. f ilt_edge_largescale parame¬ 
ter in Sections 2.1 and |2.2[ below) since the scanning speed 
of the telescope Is known and thus a length scale can be 
converted to a time scale. 


4. The AST model then identifies signihcant astronomical 
signal in the estimated map and removes its projection 
from the time series so that the noise can be measured 
accurately. 


5. Finally, the NOI model measures the noise in the residual 
time series after the removal of the AST signal. 


If the solution converges, the algorithm produces the 
hnal map. Otherwise, the process repeats itself by invert¬ 
ing the previous solution and re-estimating each model until 
convergence is achieved. 

All continuum observations presented in this paper are 
taken from the CBS survey ||Ward-Thompson et al.|2007b[) 


and were made with the PONG1800 mapping mode (Bint- 


ley et al. 20141 which produces individual circular regions 


called “pongs” with a usable diameter of ~0.5°. To map large 
regions, as in the case of Orion A South (see Section]^, cir¬ 
cular pongs are placed so their edges overlap to allow for a 
more uniform noise level in the hnal mosaic. This hnal mo¬ 
saicked dataset includes four to six repeats of every ~0.5° ob¬ 
serving held, with a higher number of repeats for data taken 
in worse weather conditions. The maps are created with 3” 
pixels and the hnal images, originally in units of picowatts 
(pW), are converted to Jy Beam“^ using the peak inten¬ 


sity c onversion factor 537 Jy pW ^ Beam ^ (Dempsey et al. 
20131 for the CBS LRl reduction and 567 JypW“ Beam“ 


(Graves et al. in prep) for the JCMT LRl reduction. For 
the CBS LRl reduction, the iterative process was termi¬ 
nated when the average pixel value changed by less than 
0.1% of the estimated map rms. In contrast, the JCMT LRl 
reduction uses a termination value of 1%. The difference 
in the stopping criteria (and other paramaters) was based 
on the specihc goals of each data reduction. A more strin¬ 
gent threshold was used for the CBS LRl because this al¬ 
lows more dihuse, extended structure to be recovered. In 
contrast, the less stringent threshold was sufficient for the 
JCMT LRl reduction as it was only important to reach con¬ 
vergence after recovering bright, compact structure which 
takes less computing time (see Sections 2.1 and |2.2| below). 
Note that the three regions presented here have a noise level 
of arms ~ 0.0038 Jy Beam“^ (for a description and analy¬ 
sis of the entire Orion A South region, see Mairs et al. in 
prep.). In the following two subsections, we discuss individ¬ 
ual choices for a subset of makemap’s configurable param¬ 
eters for each of the JCMT LRl reduction and CBS LRl 
reduction methods. 


2.1 JCMT LRl Data Reduction 


The JCMT science and computing group is currently pro¬ 
ducing for release to the wider astronomy community a 
uniform reduction and co-addition of its publicly available 


SCUBA-2 850 pm data (for more information, see Bell et al. 


(20141 and Graves et al. in prep.). This legacy release will 


consist of a) individual Hierarchical Equal Area isoLatitude 
Pixelization (HEALPix) tiles with equal area but some non¬ 
square pixels (see Gorski et al.||20d5 l reduced using a map- 
maker configuration chosen to work for a diversity of re¬ 
gions and observational types; b) these individual observa¬ 
tions coadded together to produce HEALPix tiles covering 
all the regions observed; and c) catalogues of the emission 
detected in each map. 

The makemap parameters chosen are available with the 
Starlink smurf map-maker software in a text file named 
dimmconfig_jsa_generic.lis. The parameters were devel¬ 
oped with a focus on minimising the possibility of artificial 
emission being created during the reduction process while 
still producing high quality results across a diversity of ob¬ 
servation types towards a range of astronomical regions. To 
accomplish this task, it was decided that no attempt to re¬ 
cover large-scale structure would be made. In addition, given 
the automated nature of these reductions, external masks 
rigorously dehned by an observer (like in the case of the 
CBS LRl, below) were not a viable option, and so, a more 
restrictive automasking configuration was used. The follow¬ 
ing are the most important of the makemap parameters: 


com.perarray = 1 This parameter creates a separate com¬ 
mon mode for each sub-array of SCUBA-2, which means 
that sources larger than the sub-array size (~200”) will not 
be recovered. As stated previously, accurately recovering 
compact structure is one of the main goals of the JCMT 
LRl reduction; this reduction is not suitable for an analysis 
of extended structure. 


f It. f ilt_edge_largescale = 200 This parameter Hlters 
all emission on scales above 200” and is consistent with 
com.perarray=l; it operates in the time stream and uses 
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the scanning speed to convert the length scale to a time 
scale. 

numiter = -25 & ast. skip = 5 The hrst 5 iterations are 
done without an AST model and then up to a further 20 it¬ 
erations are performed. The reduction will exit at that point 
even if it hasn’t converged. Processing of all the Orion ob¬ 
servations, however, converged. This choice of 25 iterations 
helps keep the reduction of a very large number of observa¬ 
tions to a reasonable time scale. 

ast. zero_snr = 5 Pixels which have a value of at least 5 
o'rms will be identified as astronomical signal. 

ast. zero_snrlo = 3 This parameter allows identified 
sources with pixel values of at least 5 o'rms to expand in 
area until the flux density values are 3 arms- 

maptol = 0.01 This parameter specihes when to terminate 
the mapmaking procedure. Using this reduction method, the 
process will terminate when the average pixel value in the 
map changes by less than 1% of the estimated map RMS. 


2.2 The Gould Belt Legacy Survey Legacy 
Release 1 

The GBS LRl reduction is a two-part process. In part one, a 
reduction similar in approach to the JCMT LRl reduction is 
run, i.e., flux is assigned to the AST model based on pixels 
with high signal-to-noise-ratios. This is referred to as the 
“automask” reduction. 

Part two is an additional step which is employed in 
the GBS LRl, but not the JGMT LRl reduction. After the 
automask reduction is performed, the individual maps are 
coadded for a higher SNR and the resulting image is used 
to dehne regions of likely emission. The boundaries drawn 
around the signihcant signal become the user-defined “exter¬ 
nal mask” and this mask is used to perform a second round 
of data reduction to better recover faint and extended struc¬ 
ture. In this re-reduction, instead of basing the AST model 
on pixels which achieve a specific SNR value, the pixels 
within the external mask boundaries are defined as contain¬ 
ing astronomical signal. This allows for a more well-defined 
masked area around structure which we are confident is real 
to be included in the AST model when compared using a 
single observation. 

Major differences from the GBS Internal Release 1 re¬ 
duction method (IRl; see Hatchell et al.|2013[ Buckle et al. 


2015 Rumble et al. 2015| and Pattle et al. |2015 1 include 
smaller pixel sizes (GBS LRl pixels are half as wide as GBS 
IRl pixels, allowing for better characterization of small-scale 
emission as well as more accurate peaks and positions of 
compact sources), and additional filtering of the raw data 
which better prevents the appearance of large-scale noise 
features in the reduced maps. Note that despite the at¬ 
tempts to minimize noise and reconstruct the diffuse emis¬ 
sion, there are still challenges in retrieving all the signal for 
large sources (see Section |^. The following are the most 
important parameters to be compared with the JGMT LRl 
reduction, above. Unless otherwise stated, these parameters 
are the same for both the automask reduction and the ex¬ 
ternal mask reduction: 


com.perarray = 0 No separate common mode for each sub¬ 


array of SGUBA-2 is created. Thus, sources with sizes up 
to the full array size (~400”) can be confidently recovered. 
Sources with sizes approaching and exceeding this value will 
have large-scale features subtracted by the common mode 
model. 

fIt. filt_edge_largescale = 600 This parameter, oper¬ 
ating on the time stream, filters all emission on scales above 
600”. This value is three times the size of the JGMT LRl 
reduction emission scale. 

flt.filt edge_largescale_last = 200 Only on the last 
iteration, the emission outside the AST mask is Hltered at 
200” instead of 600”. This parameter was defined to help 
suppress uncertain emission structure outside the masked 
regions. The detection of extended emission and the accu¬ 
racy of our calibration are trusted within the masked re¬ 
gions but outside the mask boundaries we cannot be con- 
hdent that the extended structure is real. By filtering the 
unmasked data harshly on the final Iteration, we remove 
the large-scale features and uncover any underlying, robust, 
small-scale sources. 

numiter = -300 ast. skip = 5 The first five iterations 
are done without an AST model, and up to a further 295 it¬ 
erations are performed. The reduction will exit at that point 
even if it hasn’t converged. All processing of the Orion ob¬ 
servations, however, converged. 

ast.zero_snr = 5 Pixels which have a value of at least 5 
O'rms will be Identified as astronomical signal. This only ap¬ 
plies to the automask reduction as the external mask reduc¬ 
tion dehnes what will be considered as astronomical signal 
using clear physical boundaries. 

ast. zero_snrlo = 0 This parameter does not allow sources 
identihed as astronomical signal to extend to lower flux den¬ 
sities. Instead, a pixel must have a value of at least 5 Orms 
to be included in the AST mask in the automask reduction. 
While the final reduction employs a user-dehned external 
mask, the boundaries of that external mask will be based 
on the sources identihed by the automask. 

maptol = 0.001 This parameter specihes when to ter¬ 
minate the mapmaking procedure. Using this reduction 
method, the process will terminate when the average pixel 
value in the map changes by less than 0.1% of the estimated 
map RMS. 


3 DATA REDUCTION COMPARISON IN 
ORION A SOUTH 

In Figures to 1^ we compare two reductions using repre¬ 
sentative areas of the GBS-dehned “Orion A South” region 
by resampling the JGMT LRl HEALPix images to match 
the GBS LRl pixel size and projection. Orion A South is a 
2.2° X 3.1° subsection of the Orion cloud complex, an active 
star-formation site approximately 450 pc from the Sun (see 
Johnstone fc BalIy|2006[|AlIen fc Davis|2008| and |Davis et al.| 


2009 


for more information). The e ntire Orion cloud co mplex 
has a mass in excess of 2 x 10® Mq (Wilson et al. 2005 \ and is 


comprised of two individual molecular clouds: Orion A and 
Orion B | Reipurth||2008 


Gomplex; 


Megeath et al. 
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2012 and references therein). The 
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northern section of the Orion A Cloud (Orion A North) 
is home to the Orion Nebula and the well-known integral 
shaped filament (Johnstone & Bally 1999| ). The Orion A 
South region is also a target of interest and, although less 
complex and dense, it exhibits several different stages of star 
formation (see Mairs et al. in prep, for more information on 
Orion A South). 

From Figures to we see that the qualitative differ¬ 
ence between the JCMT LRl and the CBS LRl reductions 
rests in the extended emission. As expected, compact struc¬ 
tures present in the data are well accentuated in the former. 
Indeed, if one’s goal is simply to determine where compact 
emission exists in the map, this style of data reduction is 
entirely reasonable. One must, however, be cautious when 
analysing the JCMT LRl maps any further as the total flux 
density present in a given region will likely be underesti¬ 
mated even in tight boundaries around the brightest sources 


(see Section 3.1.2 below). 

The CBS reconstruction gives a more accurate picture 
of the large-scale structure, as required for the science goals 
of the consortium. Beyond simply identifying where emis¬ 
sion exists, the intent here is to recover the full emission 
structure. The less drastic filtering parameter and the man¬ 
ual external masking process allow, for example, much more 
extended emission associated with a given object to be in¬ 
cluded in stability calculations, providing the opportunity 
to characterise the wide varieties of dense gas/dust mor¬ 
phologies seen in all regions of the survey. In comparison, 
the JCMT LRl reduction’s recovery of fractional extended 
emission inhibits these goals. 

In the bottom left panels of Figures to we present 
maps depicting the JCMT LRl reduction subtracted from 
the CBS LRl legacy reduction for each representative region 
of Orion A South. With only a few exceptions (blue circles 
have been drawn around prominent examples in Figures 
and|^, we see that the compact structure has been almost 
completely canceled out and only the extended emission re¬ 
mains. This situation is not unexpected since both reduction 
techniques process compact objects in the same manner. In¬ 
stead, the significant differences are found in the large-scale 
emission. 

To see the similarity of the compact structure between 
the two reduction techniques more easily, the bottom right 
panel of each Figure shows the flux density through a slice 
of each region. The solid green line shows the GBS LRl in¬ 
tensities with the JCMT LRl intensities subtracted while 
the dashed black line shows only the JCMT LRl intensities 
across the same cut. The physical positions of each slice is 
represented in the top left and bottom left panels of each 
figure by the dashed lines. Indeed, the features seen are eas¬ 
ily recognisable by comparing the image to the intensity 
plots. We see in many cases that the peaks in the JCMT 
LRl image are fainter than the peaks in the GBS LRl im¬ 
age. This is because the GBS LRl data reduction recovers 
more of the underlying flux density when using an external 
mask covering a larger area than the automatically detected 
emission in the JCMT LRl map. Evidently, the better re¬ 
covery of this extended structure occasionally increases the 
flux density of small scale features associated with the par- 


sured in the GBS LRl image. The amount in which the flux 
density increases depends on the properties of the larger- 
scale parent source. This difference can cause discrepancies 
where the JCMT LRl peak intensities are as low as 40% 
of the GBS LRl peak intensities in the observed regions of 
Orion A South. In most cases, however, the JCMT LRl peak 
intensities are ~60% of the GBS LRl peak intensities (see 
Figure]^. The brightest peaks have the smallest differences. 
Note that fainter, compact, isolated peaks also appear to 
have a bias to slightly lower intensities in the JCMT LRl 
map due to the harsh filtering parameter and less thorough 
masking procedure. This effect, however, is not as strong as 
the pedestal caused by the recovery of extended structure. 

Compact objects near the boundaries of HEALPix tiles 
can deform and elongate. The blue circle in the bottom left 
panel of Figurehighlights a peculiarity resulting from con¬ 
verting the HEALPix projection employed by the JCMT 
LRl reduction to the tangent plane projection used by the 
GBS LRl. The oblong shape of the circled source could po¬ 
tentially alter the position of a detected peak in this region. 
We note, however, that this event is uncommon and we only 
see this one example in the entire Orion A South map. 


3.1 Quantitative Differences between the GBS 
LRl and JCMT LRl reductions 

3.1.1 Structure Identification Algorithm 

To compare the two reductions fairly, we need to identify 
structure in a consistent manner across each of the two 
maps. This situation is where the GBS LRl’s diffuse, ex¬ 
tended emission presents us with a challenge. There are 
many structure identification algorithms freely available to 
apply to data such as Orion A South and each program will 
break up large-scale structures or amalgamate smaller ones 
in different ways (see|Stutzki fc Guesten[1990)|Williams et ah] 


1994[ iRosolowsky et al.|2008| [Men’shchikov et al.|2012[ and 


Berm 2015 for algorithm examples). We note that several 


individual JCMT LRl “sources” can reside within a single 
GBS LRl reduction “source” no matter which algorithm is 
used to identify structure, due to the lack of extended emis¬ 
sion in the former defining smaller areas of significant signal. 

Since we are specifically interested in the differences be¬ 
tween compact structure and extended structure separately, 
for this analysis we employ the JGMT Science Archive algo¬ 
rithm jsa_catalogue found in Starlink’s PICARD package 
I Gibb et al.|2013 (. The catalogues generated with this algo¬ 
rithm are based around the concept of islands (or extents) 
and peaks (blue contours and magenta Xs, respectively in 
Figure]^. The routine was designed to do a good job of cat¬ 
aloging all regions where emission was strongly detected (the 
islands), and then to provide a list of local peak positions 
within each island to categorise the nature of the emission. 
Source catalogues based on this approach will be released 
along with the JCMT LRl maps (Graves et al., in prep.). 

The jsa_ catalogue routine was used to identify regions 
of contiguous emission above a minimum number of pixels (9 
pixels), using an SNR cut of 5 (the noise level in these maps 
is ~ 0.0038 Jy Beam“^). These regions are found using the 


ent source, resulting in features like that which is indicated 

FellWalker algorithm (Berry 

20151 

by the blue circle in Figure The residual peak intensity 

cupid package (^Berry et al. 

2007f 


left over in this example is 12% of the peak intensity mea- 


20151 as implemented in the 


then their shapes are extracted from the data map so that 
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Figure 1. Region 1: The first of three representative regions of Orion A South. Top left: 850 pm SCUBA-2 image reduced with the 
JCMT LRl reduction parameters. Top right: 850 pm SCUBA-2 image reduced with the GBS LRl reduction parameters including the 
external mask. Bottom left: The GBS LRl map subtracted by the JGMT LRl map. The blue circle indicates a peculiarity in the map 
due to realigning the HEALPix projection of the JCMT LRl reduction to the tangent plane projection of the GBS LRl reduction. 
Bottom right: The intensities of the JGMT LRl map subtracted from those of the GBS LRl map across the positions corresponding to 
the dotted line shown in the bottom left (green) and the intensities of the JGMT LRl map across the same coordinates (black). 


the total flux density, average noise across the region, and 
the flux density at the peak can be calculated. These regions 
are identified as the islands or extents of emission. 

Within each of these detected regions, the routine then 
searches for peaks, again using the FellWalker algorithm. 
The routine is configured to identify a pixel as a peak if it 
has a) a peak value of magnitude 5 times higher than the 
average noise in the island; b) a minimum dip of magnitude 
5 times grater than the average noise towards the nearest 
larger peak; and c) more pixels than the minimum of 9 as¬ 
signed towards this peak (to avoid noise clumps). The peaks 
are identified only by their position and brightness. No at¬ 
tempt is performed to assign a shape to them. Selecting this 
set of parameters was a result of extensive testing of the Fell- 
Walker algorithm on many different SCUBA-2 maps as well 


as maps which only included artificial structure. There is no 
standard set of parameters as the algorithm can be tuned 
differently depending on the specific scientific question of in¬ 
terest. More generally, the debate about which boundaries 
to draw around different structures, especially in regions 
with complex emission, has been present in the literature 
for many years (see, for example, [Pineda et al.|[^09[ ) and 
it should be performed carefully and with an approach that 
addresses specific goals. 

3.1.2 Properties of the Structure in Each Reduced Map 

For both reductions, Tablej^shows the number of peaks and 
the number of islands detected in each region shown in Fig- 
ures[^to[^along with the centre coordinates and sizes of each 
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Figure 2. Region 2: The second of three representative regions of Orion A South. Each panel is presented in the same manner as in 
Figure ^ 


Table 1. Summary of the three regions analysed and the number of sources found within each. The RA and DEC represent the centre 
coordinates of each region. 


Region 

RA, DEC 

Size 

Number of Peaks 

Number of Islands 




GBS 

JCMT 

GBS 

JCMT 

1 

5h38m35.632s, -7:04:39.80 

0.51° X 0.51° 

47 

51 

25 

48 

2 

5h41ml6.139s, -7:58:16.84 

0.35° X 0.35° 

35 

35 

13 

32 

3 

5h39m30.365s, -7:25:23.70 

0.37° X 0.37° 

29 

44 

12 

35 


region. Table [^provides the total, maximum, and median ef¬ 
fective radius and total flux density found in each population 
of island sources as well as the peak intensity found in each 
population of sources identified by jsa_ catalogue in each re¬ 
duced dataset. The effective radius given is the radius of a 
circle with the same area as the object of interest. 

Figure shows a comparison of the source radii, to¬ 
tal flux densities, and peak intensities identified in all three 


regions. Since the islands we identify in the JCMT LRl re¬ 
duction are smaller and more fragmented than the sources 
we identify in the GBS LRl reduction, we match our JCMT 
LRl reduction sources with the associated GBS LRl reduc¬ 
tion “parent source” so that we can compare specific struc¬ 
tures in the same location of each map. We assign JCMT 
LRl sources to a GBS LRl parent based on the position 
of the geometric centre of a given source. If a JCMT LRl 
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Figure 3. Region 3: The final representative region of Orion A South discussed. Each panel is presented in the same manner as in Figure 
^ The blue circle indicates an example of a residual peak left over after the subtraction. 


Table 2. Comparison of the identified structure in the GBS LRl and JCMT LRl reductions. Three metrics are used to compare the 
GBS LRl and JCMT LRl methods in the three regions. The areas are calculated by summing the number of pixels within a given source 
identified by jsa_ catalogue’s island catalogue, the total fiux densities are the summation of the pixel values in each source’s footprint, 
and the peak intensities refer to the sources identified by jsa_ catalogue’s peak catalogue. 


Region: Metric 

Total 

Maximum 

Median 


GBS 

JCMT 

CBS 

JCMT 

GBS 

JCMT 

1: Area (arcsec^) 

3.8x10® 

6.1x10® 

2.5x10® 

7.4x10"^ 

1.4x10"^ 

7.2x10® 

1: Total Flux Density (Jy) 

142.4 

17.5 

101.8 

4.0 

0.3 

0.1 

1: Peak Intensity (Jy Beam“^) 

- 

- 

1.3 

1.0 

0.1 

0.1 

2: Area (arcsec^) 

3.3x10® 

5.2x10® 

2.4x10® 

1.1x10® 

5.5x10"^ 

6.2x10® 

2: Total Flux Density (Jy) 

149.0 

18.4 

118.4 

4.4 

2.2 

0.1 

2: Peak Intensity (Jy Beam“^) 

- 

- 

0.8 

0.6 

0.2 

0.1 

3: Area (arcsec^) 

3.1x10® 

6.3x10® 

2.7x10® 

1.4x10® 

3.3x10^ 

6.4x10® 

3: Total Flux Density (Jy) 

127.7 

26.9 

119.2 

10.0 

0.8 

0.1 

3: Peak Intensity (Jy Beam“^) 

- 

- 

2.6 

2.4 

0.2 

0.1 
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Figure 4. Islands (blue contours) and peaks (magenta Xs) identified by the FellWalker algorithm for both reductions tested using Region 
1 as an example. See text for information on the basic algorithm parameters used to identify structure. 


source centre falls within the boundaries of a GBS LRl is¬ 
land, we consider those emission structures to be associated. 
There are only a few cases where a JCMT LRl source is not 
associated with a GBS LRl source or vice versa but these 
are small, faint objects which have little bearing on our final 
results. Thus, for the radii and total flux densities, we only 
include objects in the Figure if they have a counterpart in 
each image. When comparing the peak locations between the 
GBS LRl and JGMT LRl images, the overwhelming major¬ 
ity of detected peaks lie within one pixel of their counterpart. 
Notably, there are relatively small numbers of isolated GBS 
LRl peaks and JCMT LRl peaks across the entire Orion A 
South region, this is simply due to the degree of smoothness 
of the diffuse emission recovered by each respective reduc¬ 
tion (see Section 3.1.31. As Table [^shows, region 3 has the 
largest disparity in peak number due to the manner in which 
the bright, diffuse emission is broken up by the JCMT LRl 
reduction. Many borderline peaks are identified that would 
not have been if the underlying continuous structure was 
more visible. The results presented here only include peaks 
which are associated with the same emission in each map. 


In Figure the JCMT LRl reduction source property 
is on the ordinate and the GBS LRl reduction source prop¬ 
erty is on the abscissa. For the radii and total flux densities, 
we plot each of the JCMT LRl reduction sources associated 
with a given GBS LRl source in black and we sum all of 
the associated JCMT LRl sources and plot that as a red 
plus sign. With this approach we can compare how much 
emission is being found in each region in a fair way. For 
the peak intensity, we directly compare jsa_ catalogue’s cat¬ 
alogue detections found in each map produced by the two 
data reduction methods. The blue lines in Figurej^all show 
a one to one relationship and the green lines have a slope 


of unity at the shown percentage of the GBS LRl values. 
The latter is meant to give an indication of the represen¬ 
tative difference between the two reductions based on each 
measurement. 

Since the jsa_ catalogue algorithm gives us information 
on both the peaks as well as the extended structures within 
a given map, we discuss each of these aspects individually. 


3.1.3 Peaks 

The GBS LRl and JCMT LRl reductions have similar com¬ 
pact structure in regions without extended flux (see, for ex¬ 
ample, the isolated source to the left of the main structure in 
the top two panels of Figure]^. This similarity is because the 
main differences in the data reduction parameters deal with 
structure that is much greater than one 14" beam. For exam¬ 
ple, the f It. f ilt_edge_largescale parameter (see Section 
chosen by each reduction impacts recovered emission on 
the scale of > 200". As well, the larger number of iterations 
employed in producing the GBS LRl maps will not have a 
large effect on the bright, compact structure, but that pa¬ 
rameter is set for the recovery of faint emission seen at larger 
scales. 

Since more extended emission is recovered in the GBS 
LRl reduction than in the JGMT LRl reduction, more faint 
locations are detected as individual sources. The JCMT LRl 
reduction, however, produces maps which divide up large 
structures into individual, compact components. The over¬ 
all effect is that the number of peaks detected in the GBS 
LRl and the JCMT LRl maps are very similar. In two of the 
regions, there are more JCMT LRl map peaks since loca¬ 
tions that contain diffuse emission in the GBS LRl maps will 
not meet the peak detection criteria (but do meet, in some 
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log[ GBS LR1 Peak Intensity (Jy Beam*') ] 

Figure 5. Comparing Reff (the radius of a circle with the same area as a given source identified by jsa_ catalogue’s island catalogue), 
total flux density, and peak intensity between the two data reductions in the three representative regions of Orion A South. In the 
former two, all the JCMT LRl reduction islands associated with a given GBS LRl island are plotted in black. We sum the effective 
radius and the total flux density of all of the associated JCMT LRl islands, respectively, and plot the total as a red plus sign. For the 
peak intensities, we plot the peak sources with the maximum flux density identified by jsa_ catalogue’s peak catalogue within a given 
associated island. The blue (solid) lines show a one to one relationship and the green (dotted) lines have a slope of unity at the shown 
percentage of the GBS LRl values. 
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cases, the island detection criteria) whereas the more frag¬ 
mented JCMT LRl maps will contain borderline detections 
in these regions (see Table [^. It is clear from the statistics 
presented in Tablej^that, in general, the peak intensities are 
consistently higher in the GBS LRl map in every region. As 
seen previously (recall the bottom right panel of Figure]^, 
the externally masked GBS reduction changes the amount of 
flux occasionally on the small-scale as well as the large-scale 
in regions with extended emission (see Figure]^. Evidently, 
going deeper into the GBS LRl image with an external mask 
and allowing that structure to grow raises the peaks by a 
pedestal. The pedestal is not constant, however, because it 
depends on the larger-scale, fainter emission structure. This 
behaviour is made clear once more in Figurej^which shows a 
positive correlation between the brightest JCMT LRl peak 
intensities and their associated GBS LRl peak intensities in 
log-log space. Indeed, most of the brightest JCMT LRl peak 
intensities below 1 Jy Beam~^ are only ~60% as bright as 
their GBS LRl counterparts because the diffuse background 
is not included in the former. As the peaks become brighter 
in each map, they also become more similar to one another. 
No aperture fitting is performed in this study, however. We 
perform no background subtraction for the flux extraction 
to account for large-scale structure. 

3 . 1.4 Extended Structure 

We see many similarities between the three regions shown in 
Figures to In fact, when using jsa_ catalogue to investi¬ 
gate the significant structure in the maps, it was clear that 
the identified sources in each JCMT LRl image largely trace 
the same structure as the corresponding GBS LRl image. 
The lack of diffuse emission connecting the bright, compact 
peaks in the JCMT LRl map, however, causes jsa_ catalogue 
to divide up the larger structures seen in each GBS LRl map 
into many smaller ones. Therefore, between the two reduc¬ 
tion methods, the number of “extended” sources identified 
in each region is always higher for the JCMT LRl map (see 
Table [^. The lack of extended emission in the JCMT LRl 
map is shown clearly in the relative area occupied by iden¬ 
tified sources in Orion A South. The total flux density of 
all the identified sources is closely related to the area since 
the structure identification algorithm identifies roughly the 
same peaks in both maps. 

In general, we see that the JCMT LRl image sources 
are smaller, as expected, except for the smallest (sometimes 
faint and spurious) sources present in both reduced images. 
A positive correlation is followed in radii in log-log space 
but it is not one to one. We find that the summed radii of 
all JCMT LRl islands detected within one GBS LRl island 
are typically ~45% those of the GBS LRl island radius (as 
seen in Figure]^. This corresponds to a total JCMT LRl 
island area of only ~20% of the GBS LRl island area which 
is exactly what we see in the total flux density relationship 
of Figure 


4 COMPLETENESS TESTING 

To test how well each data reduction method preserves struc¬ 
ture, we study the effect of makemap on artificial Gaussians 
with known properties. We construct and “observe” a broad 


a 


15 ' 


b 


15' 


c 


15' 


Figure 6. Top, a: An example Gaussian grid. Here, each Gaus¬ 
sian has a FWHM of 7 beams and a peak of 9 arms- The 
constructed grids are spaced accordingly for the given Gaussian 
FWHM. When the Gaussians are smaller, more sources are added 
to the noise field. Middle, 6: The field nearly devoid of structure 
in which the Gaussians were inserted. Bottom, c: The final map 
depicting the 7 beam FWHM, 9 arms peak Gaussians combined 
with the noise field using the GBS LRl reduction method. 


range of Gaussians by varying the sources’ full widths at 
half maximum (FWHM) and peak intensity values. In this 
section, we approximate 1 beam to be 15”. For each FWHM 
and peak value, we insert a grid of Gaussian sources into a 
pure noise field via makemap’s fakemap parameter (see Fig¬ 
ure]^. In this manner, the Gaussian grids are added to the 
raw time stream of the data and are subjected to the usual 
data reduction steps. We inserted 32 sets of Gaussians in to¬ 
tal with FWHMs of 1 beam (15”), 3 beams (45”), 5 beams 
(75”), and 7 beams (105”). For each of these FWHM values, 
we use a series of peak intensity values between 3 arms and 
25 ^rms • 

The field in which we insert these regular grids of Gaus- 





Data Reduction Comparison 13 


sian sources is an 850 /rm field nearly devoid of structure ob¬ 


tained by the Cosmology Legacy Survey (CLS; Geach et al. 
20131 in the PONG1800 mapping mode (see Figure [6j). This 


field is a 42 minute integration of a circular region with a 
diameter of ~0.8 degrees observed on September 30*^, 2013; 
it is a small subset of the CLS data in this region and, thus, 
it provides a representative sample of the noise in typical 
SCUBA-2 images in this configuration. We identify struc¬ 
ture in the reduced dataset by employing a Gaussian curve 
fittei|^ at the known position of the sources located away 
from the noisy edges of the map. By a simple visual anal¬ 
ysis, it appears that artificial structure outside the inserted 
Gaussians is effectively non-existent. To use the most robust 
method of returning the true structure, we provide the fidu¬ 
cial Gaussian location, peak intensity, and FWHM values 
as initial “guesses” for the peak coordinates, peak intensity, 
and the size of the expected structure to the fitting algo¬ 
rithm. We then compare the obtained output properties of 
each fitted Gaussian to the nominal input properties. 

We run each Gaussian grid run through makemap three 
times. In the first case, we use the standard GBS LRl au- 
tomasking procedure (see Section EH for tests of different 
automask parameters). Recall, however, that the GBS LRl 
reduction employs a user-defined external mask based on 
this type of automask for the final product. The second 
time, we apply an external mask which covers only half of 
the Gaussians in the image. The external mask we use is 
not generated using the normal GBS LRl external masking 
procedure. Instead, we use a simple checkerboard pattern of 
square masks, where each mask is 20 beams (5') on a side 
and centred upon every second Gaussian source (see Figure 
[^for an example and Section [5.2| for the effect of changing 
the size of the masks). For the final execution of makemap, 
we reduce the data using the JGMT LRl reduction param¬ 
eters. 

With these three sets of maps, we are able to compare 
the output properties of the GBS LRl automasked Gaus¬ 
sians, the Gaussians lying inside and outside an external 
mask, and the JCMT LRl reduction Gaussians. In addi¬ 
tion, we construct a fourth set of maps wherein the artificial 
Gaussians are spatially added directly to the reduced noise 
field. We call this latter set the “non-DR” case, as makemap 
was not used for Gaussian reconstruction. 

It is pertinent to note that these tests are all based 
on Gaussian structures. There is, however, a wide variety 
of clump morphologies, from nearly spherical dense cores 
to long, thin filaments. Gonsideration of the effects of data 
reduction on Gaussians should inform even those more com¬ 
plex sources. 


4.1 Results from the Gaussian Recovery 

The results of the Gaussian curve fitter are compared to 
the nominal Gaussian input values in Figures i to |16| In 
the following sections, we will address each of the measures 
used to compare between the output structure and the input 
objects: peak intensity, size, and total flux density. 


^ See SciPy’s generalised curve_fit function: docs.scipy. 
org/doc/scipy-0.14.0/reference/generated/scipy.optimize. 
curve_fit.html 








Figure 7. Top: The checkerboard pattern of the external mask 
for the 7 beam FWHM Gaussians. Black indicates the positive 
mask. Bottom: The final map after the GBS LRl reduction us¬ 
ing the checkerboard mask on the 7 beam FWHM, 9crrms peak 
Gaussians in the noise field. 


In each of the images, only the central nine Gaussians 
are included in the analysis to avoid the high noise locations 
on the edge of the image. This was done to mimic the ap¬ 
proach to the GBS LRl data itself as the trusted sources in 
each map lie far from the noisy edges. 

In Figures [8|to |16[ the ordinate represents the measured 
Gaussian properties divided by the nominal input properties 
at each peak intensity for each of the four FWHM values 
used. There are seven plot symbols used: 

1. The spatial addition of Gaussians onto the noise field 
(non-DR) are represented by black Xs. 

2. Gaussians included in the GBS LRl external mask (in 
the images which are produced including an external mask) 
are represented by dark blue circles. 

3. Gaussians which lie outside the GBS LRl external mask 
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(in the images which are produced including an external 
mask) are represented by light blue circles. 

4. Gausslans included in the GBS LRl automask (in the 
images which are produced without including an external 
mask) are represented by dark red squares. We define an 
object to be included in the automask if at least half of 
the total number of pixels within the respective Gaussian’s 
FWHM are found in the AST automask. 

5. Gaussians which have at least half of the total num¬ 
ber of pixels within one FWHM lying outside the GBS 
LRl automask (in the images which are produced without 
including an external mask) are represented by light red 
squares. 

6. Gaussians included in the JCMT LRl reduction au¬ 
tomask are represented by dark green triangles. Again, we 
define an object to be included in the automask if at least 
half of the total number of pixels within the respective 
Gaussian’s FWHM (in all directions) are found in the AST 
automask. 

7. The Gaussians which lie outside the JGMT LRl reduc¬ 
tion automask are represented by light green triangles. 

We performed all reduction tests using the same input 
Gaussian peak brightness values. In each Figure, however, 
the symbols have been slightly separated along the abscissa 
for clarity as they often overlap. 


4-.1.1 Peak Intensity 

In Figurej^ we see that small objects are well recovered even 
when they are relatively faint. As the Gaussian FWHMs 
are increased (top left to bottom right), the masked results 
converge only for brighter Gaussian peaks. For faint objects 
larger than a point source, however, the unmasked (light red 
squares and light blue circles for the GBS LRl, light green 
triangles for the JCMT LRl reduction) and automasked 
(dark red squares for the GBS LRl reduction, dark green 
triangles for the JCMT LRl reduction) cases display signif¬ 
icant deviations from the original Gaussian peak values in 
both reductions. As the Gaussian sizes are increased, the un¬ 
masked, reconstructed source properties become lower than 
those of the source inserted into the map. We also find that 
as a Gaussian becomes brighter, the given source is recovered 
more reliably. This behaviour is expected since the pixel-to- 
pixel variation is steeper for a cut across a Gaussian of a 
given size with a brighter peak. In contrast, shallow vari¬ 
ations tend to get filtered out. Thus, we note that larger 
Gaussians require brighter peaks for significant pixels to be 
identified and placed in the automask. 

The GBS LRl Gaussians within an external mask (dark 
blue circles) and the non-DR Gaussians (black Xs) display 
the most consistent results with the input peaks and the 
GBS LRl automask reduction (dark red squares) recovering 
the objects well, even for a 7 beam FWHM source (with 
a peak of at least 9 (Jrms- Larger objects will need to be 
even brighter to be detected). There are obvious sources 
within an external mask that are lying in a negative bowl in 
the more compact Gaussian grids (see the external masked 
Gaussians lying significantly below their nominal size in the 


top right and bottom left panels of Figure]^. This nega¬ 
tive bowl is part of the artificial structure introduced by an 
external mask that is too large for the astronomical signal 
present (see Section 5.21. As the sources begin to resemble 
the mask size, the negative bowl Is no longer apparent. For 
these and smaller sources which are sufficiently bright, the 
peak intensities are generally accurate to within 10% of the 
nominal value, which is similar to the expected calibration 
uncertainty, while larger and fainter sources are accurate 
to within 20% (with some exceptions). The variation seen 
in the small, faint non-DR Gaussians is simply due to the 
Gaussian curve fitter algorithm looking for the optimal so¬ 
lution. When it is faint, recovery of a Gaussian source will 
noticeably suffer more from noise variations, causing the al¬ 
gorithm to find the “best” peak. When the structure is small, 
any deviation from the true location will have greater im¬ 
pact on the peak intensity value and size of the Gaussian fit 
than the same deviation would have for a larger, smoother 
source with more pixels near the correct peak value in the 
same, central vicinity. 

We see also In Figure that the JGMT LRl reduction 
(triangles) reliably recovers compact structure The peaks 
found for large objects, however, are significantly underes¬ 
timated. In both the GBS LRl and JCMT LRl reductions, 
we see that the automasked sources resemble the nominal 
peak brightness more accurately as the peak brightness is 
increased before the object is defined to lie within an au¬ 
tomask. This behaviour occurs simply because the compact 
central region of the Gaussian has been identified as it is 
the brightest location of the source. Less than 50% of the 
pixels within one FWHM of the Gaussian centre, however, 
have been included in the astronomical signal (AST) mask 
at this point. 


4-1.2 Size and Total Flux Density 

In Figure we see results similar to those of the peak in¬ 
tensity for the fitting algorithm’s calculated object sizes. It 
is obvious that the GBS LRl and JCMT LRl automask 
reductions (squares and triangles) miss large structure, as 
expected. In both cases, however, they identify that at least 
some structure exists at locations of extended emission, 
assuming the peak brightness is high enough. The GBS 
LRl automask recovers nearly 50% of the pixels within one 
FWHM of a 7 beam FWHM Gaussian with a peak of 9 Urms, 
as evidenced by the increase in peak intensities in Figure]^ 
The JCMT LRl reduction requires a 15 arms peak for the 
same size, however. Once the external mask is applied to 
this structure in the second step of the GBS LRl reduction, 
the values which are measured resemble the original Gaus¬ 
sian properties. Nevertheless, occasional pedestals and bowls 
are found which can increase or decrease the size of a given 


source by up to 40% (see Section 5.2 for more information). 
As discussed above, the GBS LRl reduction is tuned to re¬ 
cover extended structure whereas the JCMT LRl reduction 
is tuned specifically to find compact emission. This means 
that the JCMT LRl reduction will always underestimate 
the large-scale structure present in an observed region. 

There are Gaussians that clearly lie on deviations in 
the noise background induced by the external masking pro¬ 
cedure (see Section |5.2| below for a discussion on external 
mask size). The noise field used is very uniform (see Fig- 
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Figure 8. Artificial source recovery comparison for different data reduction methods: peak intensities. The plot symbols have been 
separated along the abscissa for better legibility. The ordinate represents the measured output peak intensity divided by the nominal 
input peak intensity. Top left: Gaussians with a 1 beam FWHM. Top right: Gaussians with a 3 beam FWHM. Bottom left: Gaussians 
with a 5 beam FWHM. Bottom right: Gaussians with a 7 beam FWHM. 


ure however, in Figures and |10| we see the recov¬ 
ered compact Gaussians’ properties are overestimated and 
underestimated while the largest Gaussians’ properties only 
exhibit the former simply due to their placement in the map. 
The masks used here are too large, causing a pedestal effect 
on some sources and a negative bowl on others. This result 
reinforces the idea that large-scale noise features present in 
the map and the size of the external mask around a given 


source can play an important role in source recovery. Again, 
just as in the case of the recovered peak intensities, bright 
sources are found to be within 20% of their nominal value 
and many are within 10%. Looking to the non-DR data 
(black Xs), we see that there Is a fundamental difference in 
the large-scale structure when we compare with the exter¬ 
nal masked case (circles). Evidently, in some cases, makemap 
introduces additional structure in the image reconstruction 
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Figure 9. Artificial source recovery comparison for different data reduction methods: sizes. The ordinate represents the measured output 
Gaussian size divided by the nominal input size. The plotting style follows Figure]^ 


when large external masks are employed and, to conserve 
flux, it compensates by reducing the real structure present. 
This behaviour can potentially create bowls around the bor¬ 
ders of brighter sources. As the sources increase in size, the 
external mask’s pedestal effect becomes more apparent (see 
Figures]^ and 101. For an example of a bowl, see the 3 beam 
FWHM and 5 beam FWHM cases in Figurej^ The fractional 
importance of the pedestal or bowl, however, declines with 
input Gaussian peak strength. Separating the background 
sky signal from the astronomical signal using a tight bound¬ 
ary around significant structure is, therefore, very important 
so that makemap does not allow features in the noise to grow 
during subsequent iterations. As we explore below in Section 
|5.2[ an external mask should not be larger than twice the size 
of an emission region if the errors in the measured properties 
of the source are to remain within 10-20% of their nominal 
value (up to the size of the array footprint and/or charac¬ 
teristic high-pass filter scale). This gives the researcher an 


opportunity to safely mask a whole region with a generous 
boundary in case there is indeed faint extended structure 
that is not apparent from a simple automask reduction. 

It is important to reiterate that in the GBS LRl reduc¬ 
tion, the external mask size which is used is based on the pre¬ 


viously performed automask reduction (see Section 2.21. In 


this automask reduction, the common mode subtraction over 
the spatial scale of the bolometer array at each time step and 
the filter size acting on the time stream are the most im¬ 
portant considerations when determining the scales we can 
trust in the final images (see parameters com.perarray and 
f It. f ilt_edge_largescale in Section]^ also see |Chapin| 
et al. 20131. The full SGUBA-2 focal plane, including all 
four subarrays, is 400” x 400”. This nominally means that 
objects with sizes up to ~400” can be confidently recovered 
during the common mode subtraction at each time step. 
Large objects approaching this scale, however, can create 
a similar signal in a high percentage of bolometers causing 
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Figure 10. Artificial source recovery comparison for different data reduction methods: total flux densities. The ordinate represents the 
measured output Gaussian total fiux density (peak X size^) divided by the nominal input total flux density. The plotting style follows 
Figure 


them to be targeted as common mode, low frequency noise. 
Thus, large-scale modes can be subtracted from these struc¬ 
tures, diminishing their sizes and leading to uncertainties in 
drawing the external mask boundaries. 

In addition to this common mode subtraction, the fil¬ 
ter operating on the time stream was chosen to be 600”, 
approximately the diagonal size of the full 400” x 400” focal 
plane array. As noted in Section this filter uses the scan¬ 
ning speed of the JCMT combined with the provided spa¬ 
tial scale to subtract large-scale modes from each individual 
bolometer time series. Thus, on scales of ~400” to 600” and 
larger (depending on the external mask boundaries), signifi¬ 
cant structures are recovered but with diminished sizes and 
total fluxes due to the subtraction of large-scale modes (see 
Section 2 of Pattle et al.|2015 h The degeneracy between sig¬ 
nificant large-scale sources and the common mode, however, 
can cause artificial structure to arise in various parts of the 


ing external masks which are larger than necessary can fuel 
these degeneracies and create obvious pedestals and bowls. 
Note that in the case of the JCMT LRl reduction, we can¬ 
not reliably identify structures larger than 200” as this is 
the time stream filter scale as well as the subarray footprint 
(the common mode is calculated over each individual sub¬ 
array in this reduction) and no external mask reduction is 
subsequently performed. 

An interesting question for non-artificial sources is how 
much real structure we are missing at large-scales. An au¬ 
tomask will certainly pick up the brightest inner locations 
of the larger Gaussians, but, the correct size of the external 
mask to use surrounding this area is clearly debatable. An 
external mask which is too small will miss valuable structure 
but one that is too large will result in artificial structure, as 


map (see Chapin et al. 2013 section 4). Therefore, draw¬ 


mentioned above (also see Section 5.2 for a more quantita¬ 
tive analysis). Large Gaussians which are not included in 
the external mask will evidently be missed, especially faint 
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Table 3. Summary of the ast.zero_snr and ast.zero_snrlo pa¬ 
rameters tested. The ast.zero_snr parameter represents flux 
threshold for identifying astronomical signal. The ast.zero_snrlo 
parameter allows (or disallows if it is set to 0) identified sources 
with pixel values of at least the flux threshold defined by 
ast.zero_snr to expand in area until a second flux threshold is 
met. Bold font indicates the current GBS LRl automasking pa¬ 
rameters investigated in Section]^ 


ast.zero snr 

ast.zero snrlo 

Gaussian FWHM 

5 

0 

7 beams 

5 

3 

7 beams 

5 

2 

7 beams 

3 

2 

7 beams 


objects, and there is already an appreciable amount of un¬ 
recovered emission (~30%) missing from the unmasked and 
automasked Gaussians with FWHMs of only 3 beams (see 
the squares, triangles, and light blue circles). 


5 OTHER REDUCTION CONSIDERATIONS 


Table 4. Summary of the sizes of the square masks in the checker¬ 
board style external mask tests. The “size” indicated here is the 
length of the sides of the square external masks placed over every 
second Gaussian. Bold font indicates the original external mask 
size investigated in Section]^ 


Size of Square Patch in Checkerboard 

Gaussian FWHM 

4 beams 

7 beams 

12 beams 

7 beams 

20 beams 

7 beams 

36 beams 

7 beams 


5.1.2 Size and Total Flux Density 

A similar trend is seen in the recovered sizes and total flux 
densities (Figures |12| and[T^ as in the recovered peak bright¬ 
nesses, above. When the ast.zero_snrlo parameter is set to 
zero, much of the expected emission lies outside the AST 
mask for fainter objects. When the automask is extended 
from the identified 5 arms peaks down to a level of 2 arms 
or 3 arms, we see that more significant emission is recovered 
in the AST mask. As expected, the measured total flux den¬ 
sity shows very similar results to the size of the Gaussian 
structures recovered by the automask. 


5.1 Changing the Automask Parameters 

To explore how an automasked reduction would differ us¬ 
ing different thresholds for the AST model, we change the 
parameters ast.zero_snr (the threshold at which significant 
structure is identified) and the ast.zero_snrlo (the flux level 
of the surrounding pixels to which the identified significant 
peaks will be extended) using the GBS LRl automask re¬ 
duction method (see Tablej^. Since the automask reduction 
accurately recovers compact structures over a broad range 
of these two parameters and borderline island detections are 
found within the extended emission, we only perform this 
analysis for 7 beam (105”) FWHM Gaussians. Note that 
setting the ast.zero_snrlo parameter to 0 does not allow 
any extension to flux levels lower than the threshold defined 
by ast.zero_snr. 


5.1.1 Peak Intensity 

In Figure m we compare four different automask re¬ 
ductions. The reduction using the original GBS LRl au¬ 
tomask parameters are shown in the top left and the other 
three reductions were performed with the ast.zero_snr and 
ast.zero_snrlo parameters in the dimmconfig file changed 
to the values shown. Glearly, when the ast.zero_snrlo pa¬ 
rameter is not used, significant structure in faint sources is 
missed in the AST automask. When one allows identified 
structures with masked pixels (brightnesses above 5 arms) 
to grow down to a level of 2 arms or 3 arms, however, much 
more of the expected Gaussian brightness is recovered, espe¬ 
cially for the fainter input peaks. This improvement happens 
because when the constraint on the minimum flux included 
in the AST model is relaxed, more of the full Gaussian ex¬ 
tent is identified as astronomical signal earlier in the iter¬ 
ative map making process, allowing more of the significant 
emission to be extracted from the noise. 


5.2 External Mask Size 

Due to the observed pedestals and bowls induced by the size 
of a masked region around a source, we perform a further 
test of the external mask reductions where we change the 
size of each mask in the checkerboard (see Figure]^. Table 
[^outlines the different mask sizes explored, with bold font 
indicating the original external mask size used in the analysis 
presented in Section]^ Again, tests are performed only for 
the 7 beam FWHM Gaussians using the GBS LRl external 
mask data reduction and in this case we only observed a 
pedestal effect (the bowls arose from different Gaussian grid 
configurations). 


5.2.1 Peak Intensity 

In Figure |14[ we see that the application of different sizes 
of external masks can play a key role in the extracted data. 
For the small, 4 beam mask cases, much of the flux is sup¬ 
pressed since the Gaussian FWHM itself is almost double 
the mask size. Once the mask is larger than the Gaussian, 
however, we see that the expected peak intensity of the ob¬ 
ject is returned reasonably well, even in the faintest cases. 
For masks larger than 4 beams, we see the introduction of 
artificial structure (up to ~30%) and a higher uncertainty 
for the fainter sources as the pedestal effect increases for 
larger masks. Brighter peaks in general result in more accu¬ 
rate measurements, as expected, since more of the Gaussian 
emission will be more significant relative to the noise. 


5.2.2 Size and Total Flux Density 

Again, in the size and total flux density plots (Figures |15| and 
161 we see similar results to the peak intensities discussed 


above. In the smallest 4 beam mask case, we see tightly con¬ 
strained recoveries that miss the outer emission. This result 
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Figure 11. Artificial source recovery comparison for different GBS LRl automask parameters: peak intensities. The ordinate represents 
the measured output peak intensity divided by the nominal input peak intensity. Light red indicates that the object had less than 
50% of the pixels within one FWHM of the peak location detected in the AST mask, dark red indicates it had at least 50%. Top left-. 
ast.zero_snr = 5, ast.zero_snrlo = 0, the original GBS LRl automask parameters. Top right: ast.zero_snr = 5, ast.zero_snrlo = 3. 
Bottom left: ast.zero_snr = 5, ast.zero_snrlo = 2. Bottom right: ast.zero_snr = 3, ast.zero_ snrlo = 2. 


is expected because we are focusing on the brightest region of 
each Gaussian. A tighter mask that is large enough to cover 
the whole source (e.g., 12 beams in our case) produces the 
best results. It is clear that uncertainties grow rapidly for the 
faint sources. Also, in the 36 beam mask case, the brighter 
objects are also affected by artificial large-scale structure 
introduced by makemap. Ideally, the mask should just en¬ 
compass the source with little background flux at the edges 
to ensure none of the fainter extended emission is missed and 
that no substantial noise is included in the external mask. 
By using an external masking procedure based on a rigorous 
analysis for each individual region observed by the GBS, it 
is unlikely that the 850 fim maps produced for the GBS LRl 
reduction will suffer from mask areas overestimated by much 
greater than twice the size of the source. Thus, the reduc¬ 
tion performed with the 12 beam external mask masks in 


this section shows the closest example to the real GBS LRl 
data. The measured total flux density shows very similar 
behaviour to the size plots, as expected. 

6 DATA REDUCTION AND COMMON 
PHYSICAL MEASUREMENTS 

Evidently, the JCMT LRl and GBS LRl data reduction 
methods recover different amounts of extended structure. 
In this section, we qualitatively discuss the impact of using 
these different data reduction techniques on two common 
physical measurements in star-forming regions: 1. the core 
mass function (CMF) and 2. the derivation of the tempera¬ 
ture. Note that in any project, the reduction method should 
be chosen based upon the specific scientific goals one is re- 
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Figure 12. Artificial source recovery comparison for different GBS LRl automask parameters: sizes. The ordinate represents the measured 
output size divided by the nominal input size. The plotting style follows Figure 


searching with an understanding of the benefits as well as 
the drawbacks offered by that method. 


6.1 The Core Mass Function 


There have been many studies on the mass distribution of 
core populations derived from the flux due to dust in nearby 
star forming regions (for examples, see Johnstone et alT] 


2000b||Ward-Thompson et al.|2007a[|Enoch et al.|2008[ and 


Sadavoy et al.||2010 l. Although this measurement is intrin¬ 

sically dependent on the core identification algorithm used, 
specifically in the ways that boundaries are drawn around 
cores, the type of data reduction employed can also inform 
the final results. Along with the definition of the AST mask, 
the main reduction parameters under consideration when 
analysing the mass of extended objects are those which gov¬ 
ern the time stream filtering and the common mode low- 
frequency noise removal. 


In the optically thin (r << 1), isothermal limit, the 
dust emission traces the mass of a given object. A generally 
accepted result is that the more massive objects are not just 
brighter but also larger (see, for example, [Sadavoy et al.| 


2010 s Figure 12). With harsh filtering criteria like those 


employed in the JCMT LRl reduction, the most extended 
structure is not recovered and the AST mask consists of 
fairly tight boundaries around the most concentrated, com¬ 
pact regions of the map. The common mode noise is also 
estimated over each individual subarray, limiting the size of 
confidently detected structure to 200”. This may not result 
in an accurate assessment of the population of core masses as 
any larger objects will be “missing” flux. Thus, since larger 
objects correlate with larger intrinsic masses, the core mass 
function may be steepened. 


The GBS LRl is a much more robust reduction tech¬ 
nique to use for determining core mass functions when com¬ 
pared to the JCMT LRl assuming the same core identifi- 
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Figure 13. Artificial source recovery comparison for different GBS LRl automask parameters: total flux densities. The ordinate represents 
the measured output total flux density divided by the nominal input total flux density. The plotting style follows Figure [TT] 


cation algorithm. The larger size over which the common 
mode is estimated, the more relaxed filtering parameters, 
and the greater number of iterations performed work to re¬ 
cover much more significant structure than just the compact 
sources. The two-stage masking process is also highly bene¬ 
ficial for measurements like the core mass distribution. This 
ensures that two passes are taken at recovering the maxi¬ 
mum amount of structure, including a robust analysis that 
is not automated. To improve this even further, we have sug¬ 
gested that the next generation of the GBS data reduction 
be performed with different ast.zero_ snr and ast.zero_ snrlo 
parameters to potentially recover more faint structure (see 
Sections [ 5 ^ an d0. Of course, as previously discussed, large- 
scale modes will still be removed from the final map so it 
is very important to ensure that this is taken into account. 
Comparing the SCUBA-2 data with Herschel Space Obser¬ 
vatory data is one method that can be used for better under¬ 
standing the large-scale structure present in a given region. 


6.2 Derivation of the Temperature with SCUBA-2 

Another common measurement performed with SCUBA-2 
data involves the temperature maps produced by calculat¬ 
ing the 850/im/450^m intensity ratio, or against other data 
sets (such as the Herschel Space Observatory). In order to 
perform this measurement, the beams should be matched 
between the two data sets and the filtering parameters as 
well as the masks should be identical. The large-scale noise 
present in each individual map, however, will differ and 
this uncertainty becomes very important when comparing 
two maps.Preliminary work performed by GBS team mem¬ 
bers (Hatchell & Rumble, private communication; also see 
) based on the GBS LRl reduction sug- 
structures larger than 300” indeed leads 
between 850/rm and 450/rm. In order to 
perform confident measurements, that data must be post- 
processed with a tighter filter so that only the inner regions 
of these large objects have accurate associated temperatures. 


Hatchell et al.|2013 


gests that including 
to unrealistic ratios 
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Figure 14. Artificial source recovery comparison for different external mask sizes: peak intensities. The ordinate represents the measured 
output peak intensity divided by the nominal input peak intensity. Light blue indicates that the object was outside the mask, dark blue 
indicates it was inside the mask. Top left-. 4 beam masks. Top right: 12 beam masks. Bottom left: 20 beam masks (original). Bottom 
right: 36 beam masks. 


When determining the extent of the trustworthy tem¬ 
perature calculation, the profile of the structure itself must 
also be taken into consideration. It will be very difficult to 
measure accurate ratios on large, faint emission due to the 
subtraction of large-scale modes. Again, because of the harsh 
filtering parameter and the common mode removal (see Sec- 
tion |6.1[ above), the JCMT LRl reduction method would be 
appropriate for studying the temperatures of bright, com¬ 
pact objects. Calculating a temperature map for any emis¬ 
sion that lies a significant distance from peaked structure rel¬ 
ative to the sub-array footprint (< 200”) would be very un¬ 
certain. In addition, it is important to consider that the sig¬ 
nal to noise ratio is typically lower in the SCUBA-2 450/rm 
data set (compared to the 850jim data) so artificial structure 
is likely to be more prominent in these observations. 


7 CONCLUSIONS 


In this paper, we have presented and compared two meth¬ 
ods of data reduction: the JCMT LRl reduction (Graves 
et al. in prep.) and the CBS LRl reduction (Kirk et al. in 


prep.) which both employ the makemap algorithm (Chapin 
|et al.|20T^ executed with different configurable parameters. 
We used data from the Orion A South star forming region 
(Mairs et al. in prep.) to characterise the differences in the 
scale of the emission and the source morphology recovered 
by the two methods. We have measured source peak inten¬ 
sities, total flux densities, and radii across three represen¬ 
tative regions for each reduction method and compared the 
results both qualitatively and quantitatively. To further our 
Investigation, we created a series of artificial Gaussians (var¬ 
ied in size and peak brightness) and inserted them into the 
time domain of a pure noise field using each data reduction 
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Figure 15. Artificial source recovery comparison for different external mask sizes: sizes. The ordinate represents the measured output 
size divided by the nominal input size; note the change in the ordinate’s range from the figures above so the data points would be visible 
on all panels. The plotting style follows Figure [7^ 


method. We then recovered each source using a Gaussian 
fitting algorithm and compared the measured properties to 
the nominal input values to observe how much emission was 
preserved at each scale. Note that although inserting Gaus- 
sians into the time stream is an effective way to gain insight 
into the data reduction process, real structures present in the 
GBS maps are non-Gaussian and in fact most regions dis¬ 
play complex multi-scale structure. To summarise our main 
conclusions, we find: 

1. Both reductions recover the peak intensities of bright 
compact sources (sources with FWHMs of 3 beams or less 
and a peak brightness of at least 7 arms) to within 10- 
20% of the true value. The GBS LRl reduction also accu¬ 
rately recovers the peak intensities of the larger Gaussians 
whereas the peaks produced by the JCMT LRl reduction 
are diminished because of the stringent filtering param¬ 
eters. As expected, there is more uncertainty in the size 


and total flux density measurements for both reductions 
(see below). In general, we see more accurate results for 
objects which are both brighter and more compact. 

2. Although the JCMT LRl reduction only recovers com¬ 
pact emission present in a given region and the GBS LRl 
reduction recovers the extended structure, the two reduc¬ 
tion methods trace the same general areas of significant 
signal very closely. The JCMT LRl reduction’s diminished 
extended structure, however, causes island identification 
algorithms to break up large areas into many smaller ob¬ 
jects. Conversely, the GBS LRl reduction draws out these 
locations of extended emission and, thus, much more emis¬ 
sion can be recovered. These differences are illustrated in 
Figures to 

3. For faint objects larger than a point source, only a frac¬ 
tion of the true size (and, therefore, total flux density) orig- 
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Figure 16. Artificial source recovery comparison for different external mask sizes: total flux densities. The ordinate represents the 
measured output total flux density divided by the nominal input total flux density; note the change in the ordinate’s range from the 
figures above so the data points would be visible on all panels. The plotting style follows Figure [T^ 


inally present in an artificial Gaussian is recovered without 
an external masking procedure. The larger the object, the 
brighter it must be for an automask detection. In the GBS 
LRl automask data reduction, Gaussians with a FWHM 
of 7 beams must have a peak brightness of 5 Orma for the 
peak intensity to be measured to within ~20%. To measure 
the total flux density to within 20% of the nominal value, 
a peak brightness of greater than 25 arms is required. For 
a Gaussian with a FWHM of 7 beams, the JCMT LRl 
reduction will never result in jsa_ catalogue being able to 
measure a source’s peak intensity or total flux density to 
within 20% accuracy of the nominal value due to the in¬ 
herent filtering of large spatial scales. See Section |4.1| and 
Figures through |16| 

4. When identifying objects with the automask employed 
in the GBS LRl automask reduction, the most accurate 
Gaussian parameters are measured when nearby pixels are 


incorporated by extending to lower flux values. For im¬ 
proved reductions, e.g. GBS Legacy Release 2 (GBS LR2), 
we recommend setting the ast.zero_snrlo parameter to a 
value of 2. Similarly, the ast.zero_snr parameter should be 
set to 3. See Section [5.1 [ for details. One must, however, be 
cautious when extending the AST mask around significant 
peaks down to a lower flux threshold as artificial structure 
could be introduced around noise spikes. Fortunately, in 
the variety of automask tests performed, there was no sig¬ 
nificant evidence of any artificial structure outside of the 
Gaussians being included in the AST model. 

5. The external mask used in the GBS LRl reduction can 
increase or decrease the brightness of compact peaks de¬ 
pending on the surrounding region. Thus, in some cases 
where an incorrect external mask size is used, the JCMT 
LRl reduction and the GBS LRl reduction can even differ 
a little on compact scales (see Figures to [^. Artificial 
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structure in the data caused by a poorly sized external 
mask can act as pedestals and bowls, affecting the peak 
intensity, the size, and the total flux density measured for 
a given object. A mask that is too small will only highlight 
compact regions embedded within larger extended emis¬ 
sion. A mask that is too large will include noise variations 
in the AST mask. A reasonable external mask should cover 
a given source in its entirety and extend a small distance 
into the noise. It should, however, be less than a factor 
of 2 larger than a source which was reliably recovered by 
the automask reduction in order to achieve 20% accuracy 
in total flux density (see Section 4.1.2 for a discussion on 
trustworthy scales). In Section 5.2 we find that this is pos¬ 
sible even for faint, large sources. 
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